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Abstract 

A variational proof is provided of the existence and uniqueness of evo- 
lutions of regular Lagrangian systems. 



> 

in 

oo ! Introduction 
O 

Let Q be a smooth, finite dimensional manifold and L : TQ — ► K be a smooth 
■ Lagrangian. Evolutions of the Lagrangian system defined by L are by definition 

f- > ) the C 1 curves q: [0, /i] — + K which are critical points of the action 

S h = / Loq'(t)dt, 

J: 







subject to the constraint that q(0) and q(/i) are constant. A typical route to 
existence and uniqueness (given that L is regular) of the Lagrangian evolutions, 
is to to show that derivatives q'{t) of evolution curves q(t) are integral curves 
of the Lagrangian vector field Xe, constructed either using the Euler-Lagrange 
equations in charts, or using the Lagrange two- form u>l, the energy E, and the 
equation 

ix B UL = dE. 

In any case, standard ODE theory provides existence and uniqueness of the 
initial value problem q(0) = qo,q'(0) = vq. For a self-contained exposition, 
see p. 

Given two nearby q\, qi £ Q, does there exist a unique evolution curve q(t) 
such that q(0) — qi and q(h) — q-p. This is the local boundary value problem 
of Lagrangian mechanics. The problem crops up in a variety of situations. For 
example: 

1. If Q is a Riemannian manifold and L(v) = hg(v,v), then the Lagrangian 
evolution curves are constant speed reparameterizations of the geodesies, 
and the local boundary value problem becomes that of locating the unique 
local geodesic connecting two sufficiently nearby points. 
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2. A solution to the local boundary value problem is required to construct 
type 1 generating functions St(q%, q\) for the Hamiltonian flow, which are 
defined by 

St(q2,qi) = / Loq(t)dt 
Jo 

where q(t) is the evolution curve with q(0) = q\ 1 q(t) = qi- 

After constructing the Lagrangian flow F* E , the solution to the local bound- 
ary value problem is obtained by solving the equations 

T Q F t XE Ki ) = 92 

for v qi £ T qi Q as a function of q\,q2,t, where tq : TQ — > Q is the canonical 
projection. This may appear to be a straightforward application of the implicit 
function theorem near t = 0, q\ = q%, but that is not quite so, because the 
equation fails to be appropriately regular there. With some care, however, the 
local boundary value problem can be solved by this route 

But, first solving the initial value problem seems like a rather circuitous 
route to the solution of the local boundary value problem, especially consid- 
ering that the boundary values q\ and qi actually occur as the constraints in 
the original variational formulation for the evolution. Regular constrained opti- 
mization problems have critical points which persist as smooth functions of the 
constraint values. Why not solve the boundary value problem directly, avoiding 
an excursion into the initial value problem via ODE theory? 

The obstruction to simply getting on with the job is a fundamental one: the 
problem of finding the critical points of Sh subject to the constraint q(0) = qo, 
q{t) = q\, is nonregular at h = 0, which precisely where one wants to perturb 
from. Indeed, the objective function Sh is actually zero when h = 0. More 
seriously, the constraint 

q(t)^(q(0),q(h)) 

maps into the diagonal of Q x Q at h = 0, and hence cannot be, even formally, 
a submersion. One cannot perturb from such a degenerate landscape. 

In this article, I provide a direct variational proof of the local boundary 
value problem, using a regularization procedure, which is adapted from the 
one used in [3] for a similar problem in the context of discrete Lagrangian 
systems. The regularization procedure culminates in the replacement of the 
variational problem with an equivalent, regular one, after which readily available 
techniques used to prove the infinite dimensional Morse Lemma, and the (infinite 
dimensional) implicit function theorem, give the result. 

1 Regularization 

Assume that L: Q — > R is C r with r > 2. Since the aim is to provide a local 
result, also assume that Q is an open subset of R n , and that q G Q is given. 
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One seeks a pcrturbativc approach from t = 0, q\ = q 2 = q, so that q± and q 2 
will be near q. The regularization procedure is, step-by-step, as follows: 

1. Transform the variational problem from one for curves in Q to one for curves 
in TQ, with an additional first order constraint. The transformation is simply 
to seek critical points (q(t), v(t)) G TQ of the objective 

f h 

S h = / L(q(t),v(t)) dt, 
Jo 

subject to the constraints 

v(t) = ^, q(0) = qi, q(h) = q 2 . 

There are additional freedoms inherent in the use of curves in TQ rather than 
curves in Q, and these will be important for the regularization. 

2. Reparameterize, so that the solutions curves, which are defined on [0,h], do 
not disappear as h — > + . The curves (q(t), v(t)) , t G [0, h] are replaced by the 
curves (Q(t),V(t)), t G [0, 1] through 

Q(u)=q(hu), V(u) = v(hu), u e [0, 1]. 

For h > 0, the new curve (Q(u), V(u)) satisfies an equivalent variational prob- 
lem, which can be worked out as follows. First, substitute u — t/h and divide 
by h to obtain 

\ [ L(q(t),v(t))dt= [ L(Q(u),V(u))du. 
" Jo Jo 

Since h is constant for the variational principle, one can use the right-hand- 
side as an objective for (Q(u), V(u)) . The first order constraint transforms as 
follows: 



dt 

Multiplying this by h, the reparameterized variational principle is for curves 
(Q(u), V(u)) with values in TQ which are critical points of the objective 



S= [ L(Q{u),V(u))du, 
Jo 



subject to the constraints 

A.(Q(u))-hV(u) = 0, Q(0) = gi, Q(l) = q 2 . 

Notice that both the objective and the constraints are smooth through h = 
and that, at h — 0, the first order constraint is equivalent to the constraint that 
the curve (Q(u), V(u)) lies in a fiber of TQ i.e. the curve is vertical. 
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3. Restrict the remaining (boundary) constraints to the submanifold defined by 
the regularized first order constraint, and regularize the result. The first order 
constraint may be solved smoothly through h — by integration: 

Q(u) = Q(0) + h / V(s)ds. (1) 
Jo 

Thus the set of curves (Q(u),V(u)) may be replaced by the set of curves 
{V(it)}. Also, <72 may be replaced by q\ + hz, z e R", and then, since un- 
der (P), 

Q(l) = Q(0) + h { V{s)ds, 



the objective and the boundary constraint become, respectively, 

S h = [ ifqi+h [ V(s)ds,V(u)) du, [ V{u)du = z, (2) 



Jo \ Jo J Jo 

while h and q\ appear as parameters. 

This completes the procedure, since the variational problem (0 is formally 
regular through h = 0. Indeed, at h = 0, @ becomes finding the curves V(u) 
which are critical points of the constrained problem 

So— L(q\, V(u)) du, / V{u)du — z. 
Jo Jo 

Using a Lagrange multiplier A, the solutions are given by setting 
( qi ,V(u))SV(u)du= / \-5V(u)du 



o dv j 
for all 5V(u) i.e. 

r)T 

If L is a regular Lagrangian then this implies V(u) is constant, and the constraint 
then implies V{u) — z. Thus the solution to © at constraint value z is the 
constant curve z, and there is exactly one critical point for each constraint value. 



2 Implicit function theorem solution to the reg- 
ularized variational problem 

Consider St from 101 on the Banach space C fe ([0, 1],K") of curves V(u), where 
< k < r. Since L is C r , the Omega Lemma ( 2 , page 102) implies that 
the integrand of Sh is C r as a map into C°([0, 1], R). Since integration on 
C°([0, 1],R) is bounded linear (and therefore C°°), it follows that S h is C r , 
irrespective of the value of k. 
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I specialize the method for proving the infinite dimensional Morse lemma 0J 
[7||S] to the constrained variational problem (JTJ. One first calculates the gradient 
of S from the derivative dSh using the C 2 ([0, 1],R") weak inner product 

((V, W))= ( V ■ w. 
Jo 

That computation is as follows (for short, below Q{u) means the right side 
of O): 



dS h ([V(u)])SV(u) 
d 
de 



= -j- / L( qi + h I V(s) + eSV(s) ds, V(u) + e6V(u) ) du 



1 /r)7" i9f 
' -(Q(u),V(u))h SV(s)ds + —(Q(u),V(u))SV(u))du 



Bv 

1 BT f 1 81 

h—(Q(u),V(u))SV{s)duds + J —(Q(u),V(u))SV(u)du 

1 BT f 1 BT 

h—(Q(s),V{s))SV(u)dsdu+ / — (Q(u),V(u))SV(u)du 
oq J ov 

1 BT BT \ 

h—(Q(s),V{s)) ds + — (Q(u),V(u))\ 5V{u) du, 



from which 



BL f 1 BT f u 

VS h = —(Q(u),V(u))+h J —(Q( s ),V(s)) ds, Q(u)= qi + hj V(s)ds. 

By the same reasoning as was used to find the differentiability of Sh , the gradient 
VSh (best thought of as a vector field) is a C"" -1 map from C k ([0, 1],R") to 
C°([0, 1], R"), also irrespective of the value of k, < k < r — 1. 

The constraint of (0) (i.e. the second equation) is C°° because it is bounded 
linear, and its derivative is 

SV(u) i — > / 8V{u)du. 



The kernel of this derivative, say E , is the tangent space to the constraint 
set, and it splits C k ([0, l],R n ) orthogonally with respect to the metric ((, )) (the 
complement is the subspace of constant functions) by 

SV(u) = (sV(u) - SV{u)^j © J 5V{u) 

Sh has a critical point on the level sets of the constraint if and only if the 
orthogonal projection Pe of VSh to the kernel Eo is zero i.e. for the solutions 
to the constrained variational problem J5J), one solves 

PE V5ft(F ®Fi) = 
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for Vq near 



V =0, Vi=0, z = 0, qi 



h = 0. 



To use the implicit function theorem, one requires that the appropriate partial 
derivative of Pe VSh is a linear isomorphism. Remembering to set h = 0, that 
derivative is 



d 
1c 



dL d 2 L 



dv 2 



d 2 L 
dv 2 



(q,0)8V (u). 



If L is regular, this is a linear isomorphism of Eo, with inverse 



SV (u) 



Thus, the implicit function theorem provides neighborhoods W\ C K" x IR ra x 
M = {{ qi ,z,h)} containing (£,0,0) and W 2 C C* fc ([0, 1],R") of the constant 
map u \— * 0, and a C r_1 map ^: Wi — > W2 such that for all (qi,z,h) e Wi, 
■0(<7i , z, ft) £ C fc ([0, l],R n ) is the unique critical point in W2 of the constrained 
variational problem (J2J). By setting k — and then fc = r — 1, one can arrange 
that W2 is a C r_1 neighborhood, ^ has values in W%, and hence in the C"* -1 
curves, t/» is C r ~ 1 with the C" 1 topology, but that ip provides the unique 
solution among the C° curves in a C° open neighborhood, say 

{V{u) : \V(u)\<e}, 

of the constant curve 0. 

Now reverse the regularization. Pick an h > such that (q, 0, h) G Wx, set 

VKi = { ( qi ,q 2 ) : (ft, (g 2 - qi)/h, h)eW x } 

and define 



V>(<ri, 9 2)(*) = ft + 



■0 ( 3l>r(?2 ~qi),h ) (u)du. 



Then (g, g) 6 Wi, $( qi ,q 2 ){t) is defined for (51,92) G and i 6 [0,/i], and 
i>r qi> q 2 ){t) is a first order curve in TQ which has base integral curve a Lagrangian 
evolution. This evolution is unique among the continuous curves corresponding 
to |V(u)| < e i.e. among C 1 curves q(t) such that \q'(t)\ < e/h, so C 1 curves 
q(t) in some C 1 neighborhood of the constant curve q. 
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3 Remarks 



The regularization can be formulated in invariant terms on the manifold Q, using 
a tubular neighborhood of the antisymmetric normal bundle of the diagonal of 
Q x Q to accomplish the subtraction qz — q\. Replacing (Q(u), V(u)j with its 
TQ version V(u), the regularized variational problem, at h — 0, becomes 

Sq = LoV, / V(u) — z, tqV = constant. 
Jo Jo 

It is a pretty result that the variational principle on Q regularizes to this trivial 
one on the fibers of TQ. The map i/j is defined only for small z = ((72 — Qi)/h, 
and since ip{qi, 52) is a solution which goes from qi to (72 in time h, the velocity 
of this solution is also, approximately, (92 — Thus regularizing only at 

z = provides evolutions which correspond only to velocities near zero. This 
is unacceptable for the objective of solving the initial value problem by first 
solving the local boundary value problem, because it assigns evolutions only 
to those initial data corresponding to velocities near zero, while it is known of 
course from ODE theory that there is a unique integral curve of the Lagrangian 
vector field corresponding to any velocity. However, minor extensions of the 
above show that the variational principle actually regularizes at all z. The local 
solutions so obtained along the entire solution (i.e. z maps to the constant 
curve V(u) = z) of the regularized variational principle at h = 0, may be glued 
together using a technique the can be found for example in 0, page 97. This 
provides solutions starting at any velocity. For the same reason, this also an 
important step in the discrete Lagrangian context [3], where the discrete initial 
value problem is addressed by first solving the local boundary value problem. 
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